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ABSTRACT 

Context. In recent years, many very interesting observations have appeared concerning the positions of Wolf-Rayet (WR) stars in the 
Hertzsprung-Russell diagram (HRD), the number ratios of WR stars, the nature of type Ibc supernova (SN) progenitors, long and soft 
gamma ray bursts (LGRB), and the frequency of these various types of explosive events. These observations represent key constraints 
on massive star evolution. 

Aims. We study, in the framework of the single-star evolutionary scenario, how rotation modifies the evolution of a given initial mass 
star towards the WR phase and how it impacts the rates of type Ibc SNe. We also discuss the initial conditions required to obtain 
collapsars and LGRB. 

Methods. We used a recent grid of stellar models computed with and without rotation to make predictions concerning the WR 
populations and the frequency of different types of core-collapse SNe. Current rotating models were checked to provide good fits 
to the following features: solar luminosity and radius at the solar age, main-sequence width, red-giant and red-supergiant (RSG) 
positions in the HRD, surface abundances, and rotational velocities. 

Results. Rotating stellar models predict that about half of the observed WR stars and at least half of the type Ibc SNe may be produced 
through the single-star evolution channel. Rotation increases the duration of the WNL and WNC phases, while reducing those of the 
WNE and WC phases, as was already shown in previous works. Rotation increases the frequency of type Ic SNe. The upper mass limit 
for type II-P SNe is ~ 19.0 Mq for the non rotating models and ~ 16.8 Mq for the rotating ones. Both values agree with observations. 
Moreover, present rotating models provide a very good fit to the progenitor of SN 2008ax. We discuss future directions of research 
for further improving the agreement between the models and the observations. We conclude that the mass-loss rates in the WNL and 
RSG phases are probably underestimated at present. We show that up to an initial mass of 40 Mq, a surface magnetic field inferior to 
about 200 G may be sufficient to produce some braking. Much lower values are needed at the red supergiant stage. We suggest that 
the presence/absence of any magnetic braking effect may play a key role in questions regarding rotation rates of young pulsars and 
the evolution leading to LGRBs. 
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1 . Introduction 20 million stars is a WR star in the Galaxy. To date, 476 galactic 

„,,__, ,,, rl ^ , . , , . WR stars have be en detected, i.e. abou t 7-8% of the total Milky 

Wolf-Rayet (WR) stars are characterised by four importan t w popu i ation (fMauerhanet alhoilh . Although these stars are 

IcSher 2007) reS reC6nt ^"^ ^^^ ^ ""' ^ ^ important in astro P h y sics for man y reasons: 

1 . They are associated with young massive star regions. " The y represent an evolved state of the most massive stars 

2. They present broad and strong emission lines. (see e .g. |Schnurr et al.|2008l |Rauw et aLEJOS |Lamers et al.| 

3. They are hot (log(r e ff/K) > 4) and luminous stars | 1991| ). 

(log(L/L ( T 1 ) > 5 0) - They allow us to check the nuclear reaction chains duri ng the 

4. The chemical composition of their surface shows signs of H " and He-burning phases (see e.g.|Dessart etap|2000|). 
H-burning (WN-type) and/or He-burning processes (WNC, " The y suffer stron g mass loss > making them very inter- 
WC WO-type) esting laboratories for studying the physics of stellar 

winds. Many of them are surrounded by extended nebulae 
These features can be satisfactorily explained if WR are mas- dStock & Barlowll2010b 
sive evolved stars whose surface composition has been changed _ They ^ candidate p r0 g enitors for type Ibc supernovae 
by mass loss and/or internal mixing. The mass loss can be due (SNe; Q which give bhth t0 either a neutron star (NS) or a 
either to stellar winds and/or to the loss of the envelope through black hole (BH ISmarttll2009l) 
a Roche lobe overflow (RLOF) in a close-binary system. ' 



The estimated total number of W R stars in the Milky Way is ' In this paper, we call the sample composed of type lb and type Ic 
about 6000-6500 dShara et al.l2009l) . It means that about 1 out of SNe type Ibc SNe. 
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- They may be the progenitors of long and s oft gamma-ray 
bursts (LGRBs), or at least of part of them dWooslevll 19931 
l201lh . 

- They are important contributors to the injection of new syn- 
thesiser! nuclear species in the interstellar medium through 
their winds and possibly also through their SN ejecta, mak- 
ing them an important agent of chemical evolution in galax- 
ies. In particul ar, they may be significant sourc es of 12 C 
dMaederi 119921) . ™F dMevnet & Amoukj fc OOO) and 26 A1 
(Dearborn & Blakel [l984: Pal acios et alJl2005bh . to cite only 
a few elements. A WR star proba bly injected the 26 A1 
present at the birth of the solar svstem JjTatischeff et aL 2010; 
iMontmerle et alJl20ioHGounellell201 li 

- The presence of WR stars can be detected through the analy- 
sis of the integr ated light spectrum of very distant galaxies 
dAllenetal.lll976t IContil Il99lh . Indeed, their broad emis- 
sion lines can be observed superposed to the galactic con- 
tinuum in young active star-forming systems, making these 
stars very useful tracers of distant starburst regions. 

Amongst the questions concerning the WR stars, which are a 
source of lively debate these days, we cite the following two: 

- In what way are the single-star and close-binary chan- 
nels significant for explaining the observed WR populations 



dPodsiadlowski et al.ll992tllzzard et al.l2004tlEldridge et alj 
2008)? Do these two scenarios have different effects at var- 
ious metallicities? Depending on the answers to these ques- 
tions, the age associated with a WR population might be 
quite different, the range of initial masses evolving into the 
WR phase being likely different in the single and binary sce- 
nario. 

- What kind of SNe do the WR stars produce? If they col- 
lapse into a BH, is the latter associated with a faint or does 
it prevent any SN event? Do they give birth to a type lb 
or type Ic SN? Do WR stars produced through the single 
and binary scenario have the same SN output? Answers to 
these questions are important for the contribution of these 
stars to the chemical evolution of the galaxies, for interpret- 
ing the observed frequencies of type Ibc SNe, and to check 
the possibility that these stars are in some circumstances as- 
sociated to LGRB events, since in a handful of cases the 
spectrum of a ty pe Ic SN has been observed in association 
with an LGRB (faalamaet alj[l99i IChornockeTaTll2010t 
iBergeret al.ll201l[ and references therein). 

In the present work, we aim to make progress towards pro- 
viding answers to the above questions, analysing the results we 
have obtained in our most recent solar-metallicity grid of stellar 
models (Ekstrom et al. 2012, hereafter Paper I) . These models, 
which include the effects of rotation, are able to reproduce many 
observed characteristics: 

- the characteristics of the Sun at its present age, 

- the observed width of the main-sequence (MS) band in the 
Hertzsprung-Russell diagram (HRD), 

- the positions of red giants and red supergiants (RSG) in the 
HRD, 

- the observed surface composition changes in B-type dwarfs 
and supergiants, 

- the observed averaged rotational surface velocities. 

We aim to check whether these models are also able to reproduce 
the observed ratios of WR to O-type stars, of WN to WC stars, 
and of WNC stars (a transition stage between the WN and WC 
stages where both H- and He-burning products are observed) to 



WR stars at solar metallicity. This study will fuel the discussion 
of how rotation modifies the evolutionary scenarios in the upper 
part of the HRD. We also study the nature of the SNe arising 
from stars more massive than about 8 Mq. 

The paper is structured as follows: Sect. [2]briefly recalls the 
main physical ingredients of the models. The WR stellar mod- 
els are presented in Sect. [3] Comparisons with the observations 
are performed in Sect. |4] The nature of the progenitors of type 
Ibc SNe is the subject of Sect. [5] while Sect. [6] focuses on the 
questions of the rotation rate of pulsars and the conditions for 
obtaining collapsars and LGRB. A synthesis of the main results 
is presented in Sect. [7] 

2. Physical ingredients of the models 

The phy sical ingredients of the models are described in detail in 
iPaper li We recall the main features here: 

- The initial abundances were set to X — 0.720, Y = 0.266 
and Z = 0.014. T he mixture of heavy e lements is taken to 
be the same as in [Asplund et al. (200 5J) except fo r the Ne 
abundance, which we took from lCunha et al.l (2006). 

- The opacities and nuclear reaction rates were updated (see 
details in Paper I). 

- The convective core is extended with an overshoot parameter 
d ver/Hp = 0.10 starting from the Schwarzschild limit. 

- The outer convective zone is treated according to mixing 
length theory, with qtmlt = 1.6. The most luminous mod- 
els (M > 40 Mq) can have a density inversion near the 
surface due to supra-Eddington luminosity layers. In those 
models, the mixing length is based on the density scale: 
ttMLT = f/H p - ({a - 5V)/ Hp - 1. The turbulen ce pres- 
sure and the energy turbulent flux are included (see iMaederl 
2009, Sect. 5.5). 

- Models with rotation account for the effects of the strong 
horizontal turbulence, the vertical shear, and the meridional 
circulation as explained in IPaper 1 No magnetic field is as- 
sumed. 

- A treatment allowing the precise conservation of angular mo- 
mentum in rotating models was implemented. 

- The r adiative mass-loss rate adopted is from IVink et all 

(2001). In the domains not covered by this prescription, 

i 1 j 1 

we used that of de Jager et al. ( 1988). For red (super)giants, 

we used the R eimersl ([1975, 1977 ) formula fo r stars up to 
12 Mq, with 77 = 0.5, and that of Ide Jager et al.l (1 1988b from 
15 Mq and up for stars with log(r e ff/K) > 3. 7, or a linear fit 
on the data from Sylvester et al. (1998) and Ivan Loon et all 
(fl999l) (seelCrowther 20PJ) f or log(7W/K) < 3 . 7. The WR 
stars are compu ted with the iNugi s & Lamers (2000) pre- 
scription, or the iGrafener & Hamannl (120081) recipe in the 
restricted validity domain of this prescription. Note that it 
may occur that the WR m ass-loss rate is low er than the one 
that would result from the Vink et al. (2001 ) prescription. In 
that case we keep to the Vink et al. (2001) prescription until 
the WR prescription becomes high er. These mass-loss rates 
account for some clumping effects (Mu iires et al.ll201 H) and 
are a factor of 2-3 lower than the rates used in the standard 
1992 grid. 

- When, for massive stars (> 15 Mq) in the RSG phase, some 
points in the most external layers of the stellar envelope ex- 
ceed the Eddington luminosity of the star (Lgdd = 4-ncGM/K, 
with k being the opacity), we artificially increase the mass- 
loss rate of the star (computed according to the prescrip- 
tion described above) by a factor of 3 so that the time- 
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averaged mass-loss rate during the RSG phase for 20 and 
25 Mq stars is o f the order of the mass-loss rates estimated 
bv lvan Loon et al.1 (120051) for RSGs. 



3. The WR stellar models 

3.1. HR diagram 

We used the following criteria to determine the type of the star at 
a given time dMevnet & Maederll2003HSmith & Maederill99ll) : 

- A star with log(r e ff/K) > 4.0 and a surface hydrogen mass 
fraction Xu < 0.3 is considered as a WR statQ 

- A star with log(r e ff/K) > 4.5 that is not a WR star is an 
O-type star. 

- A WR star with a surface hydrogen mass fraction Xu > 10" 5 
is a WNL star. 

- A WR star without hydrogen and with a carbon surface abun- 
dance lower than the nitrogen abundance is a WNE star. 

- A WR star without hydrogen and with a carbon surface abun- 
dance higher than the nitrogen abundance is a WC or a WO 
star. If the surface ratio ^^ (in number) is less than 1, the 
star is classified as a WC; otherwise, it is classified as a WO. 

- A WR star without hydrogen, and with a surface ratio -S 
between 0.1 and 10 is a WNC star. 

The O-type, WNL, WNE, WC, and WO phases are exclusive. 
The WNC phase lies between the WNE and WC phases. We 
point out here that the observational classification of WR stars is 
based on spectroscopic characteristics and not on surface abun- 
dances. A firm classification of our models would therefore re- 
quire computing the output spectrum. While some attempts in 
that direction have alread y been made (Schaereret al. 1996a b; 
Scha erer & de Koterll 19971) . this is by far not yet a standard pro- 
cedure. We decided here to follow the usual procedure most 
often found in the literature to connect WR spectral types to 
surface abundances. Below we present some numerical exam- 
ples that illustrate how numbers change when different rules are 
adopted for the connection between surface abundances and WR 
subtypes. In Fig.Q]we show where each WR phase occurs in the 
HRD for the models between 20 and 120 Mq. 

Comparing the rotating and non-rotating models, we note the 
following differences: 

- The MS rotating tracks for the most massive stars (M > 
60 Mq and Vj n ; > 350 km s _1 ) follow a nearly homogeneous 
evolution during the MS phase (Maeder 1987). As a numeri- 
cal example, the mass fraction of hydrogen at the stellar sur- 
face is 0.32 in the rotating 120 Mq model when the mass 
fraction of hydrogen at the centre is 0.28. We note that al- 
though these models remain quite compact and hence blue 
during the MS, their surface velocity decreases a lot. Indeed, 
at the end of the core H-burning phase, the surface velocity 
is only a few km s -1 for stars with initial masses equal or su- 
perior to 60 Mq. This is caused by the high mass-loss rates 
undergone by these stars. 

- Without rotation, only the 120 Mq enters the WR phase at 
the end of the MS, while with rotation, the 60 Mq and more 
massive stars already enter into the WR phase during the MS 
phase. 



2 The duration of the corresponding WR phase is not very sensitive to 
th e choice of this limi t if it r emains in the range 0.3 - 0.4 as mentioned 



by Meynet & Maedeil J2003h . 



- A small range of initial masses go through an RSG 
stage before entering the WR regime. These models reach 
log(r e ff/K) ~ 3.5 at the RSG stage, and then, due to the 
strong mass loss, evolve back to the blue part of the HRD. 
They finally become WR stars at the very end of their evolu- 
tion, ending as a WNL or WNE stars with an effective tem- 
perature (at the border of the hydrostatic core) of roughly 
log(r e ff/K) ~ 4.5. Rotation restricts the mass domain of the 
family of tracks evolving into a WR phase after a RSG phase 
(see Fig. 0. As a consequence, rotation decreases the max- 
imum luminosity of the RSGs from 5.7 without rotation to 
about 5.4 with rotation. The last value agree s better with the 
upper limit of RSGs, which was found by iLevesque et al.l 
(2005) to be between 5.2 and 5.4. 

- The minimum luminosity reached by the WN and WC stars 
is given by the minimum initial masses of the stars that go 
through WNL and WC phases (see Table 2). They are low- 
ered by rotation: for non-rotating models, the minimum lu- 
minosity of WNL and WC stars is 5.37 and 5.5, respectively, 
while for the rotating models it is equal to 5.25 and 5.35, re- 
spectively, i.e. in the same range as the maximum luminosity 
for RSG stars. 

- Both rotating and non-rotating models predict that the least 
luminous WR stars are WNL stars. The least luminous stars 
are produced by the least massive star that just succeeds to 
enter the WNL phase at the end of its evolution. Note that if 
the mass-loss rates are allowed to be higher than accounted 
for in the present grid, there might be situations where the 
least luminous stars would be WC stars originating from 
stars well above the minimum initial mass of stars evolv- 
ing into a WR stage. This is what was obtained in the grid 
with enhanced mass loss by Maeder & Meynet (1994), for 
instance. The nature of the less luminous WR stars therefore 
depends quite closely on the mass-loss rate history. 



3.2. WR lifetimes and mass limits 

The WR lifetimes are indicated in Table [TJ and are plotted as 
a function of the initial mass in Fig. [2] The minimum masses 
leading to a WR star (or a subtype) M '™ are given in T able [2] 
The mass limits are determined as in Georg y et al.l (120091) . 

Qualitatively, t he results do not differ from th ose presented 
in previous papers (Meynet &Maederl l2003. 2005). The general 
trend is an increase of the WR phase lifetime with increasing 
mass and rotation. The most spectacular increase concerns the 
duration of the WNL phase for rotating models, which is several 
times longer than for non-rotating ones. On the other hand, the 
WNE phase almost disappears for rotating models. Again, this 
is a consequence of rotational mixing. To understand this, we 
have to keep in mind a few specific points. First, a WNE phase 
implies some pure He layers, which obviously appear only when 
H-buming is terminated in those layers. Second, the width in 
mass of the pure He-rich region, which when uncovered gives 
the WNE phase, is decreasing rapidly with time because of the 
fast growing He-burning core that enriches the central layers in 
carbon and oxygen. Third, mixing prolongs the WNL stage, as 
just seen above, and consequently leaves more time for the star 
to evolve into the core He-burning phase, which in turn reduces 
the mass of pure He-layers. This makes the WNE phase very 
short in rotating models, or even more generally in any models 
with some efficient internal mixing. 
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Table 1. O-type star, RSG, and WR lifetimes for the most massive models. 
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Notes. The initial masses are given in Mq. Lifetimes are given in years. 
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Fig. 1. HRD of the massive models from 20 to 120 Mq with the different types/phases marked in colours (O-type: blue; neither 
O-type nor WR: red; WNL: green; WNE: purple; WC: cyan). Left: Non-rotating models. Right: Rotating models. We plotted the 
effective temperature at the surface of the hydrostatic core. The endpoints of the tracks are indicated by a circle. 



4. WR star populations: comparisons with the 
observations 



to their binary status, these two stars are not plotted in Fig. [3] 
Looking at Fig. [3] a few points deserve to be emphasised: 



4.1. Positions in the HR diagram 



Observed positions of some galactic WN (lHamann et al.| [2006) 
and WC ( S an der et al.l 120 1 2) stars are indicated in Fig. [3] The 
luminosity of the most luminous WNL star in the sample of 
Hamannet al. (2006) (not shown in Fig. [3] se e below) would 
requir e an initial mass of the order of 300 Mq. ICrowther et al.l 
(2010) have recently claimed to have observed such massive ob- 
jects. However, WR 25 (or HD 93162), the most luminous star, 
is a binary star with a spectroscopically determined orbital pe- 
riod of 207.8 d(Raassen_etalj2003). T he secon d-most luminous 
star, WR 22, is an eclipsing binary. iRauw et al. (1996) estimated 
a lower limit for the actual mass for this star to be 72 Mq. Owing 



- The WNL, WNE, and WC present some kind of decreasing 
luminosity sequence. There are of course overlaps between 
the different WR types; however, just by eye one can notice 
immediately from Fig.[3]that the averaged luminosity of each 
class decreases from the WNL to the WC stars. The WR stars 
with the highest luminosity are therefore the WNL stars (the 
three upper points are also those with the highest abundances 
in hydrogen, with values between 0.35 and 0.53), and the star 
with the lowest luminosity is a WC star. 

- The range of effective temperatures (taken here at the sur- 
face of the hydrostatic core, i.e. not accounting for the optical 
thickness of the wind) of the WNL stars is relatively narrow 
(between 4.6 and 4.8) above a luminosity equal to about 5.8. 
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Fig. 2. Lifetimes in the RSG phase (defined as stars with log(r e ff/K) < 3.66, see Eldridge et al. 2008) and in the different phases of 
WR stars. Left: Non-rotating models. Right: Rotating models. 



Below, it widens and extends from 4.6 up to slightly more 
than 5. 
- The ranges of effective temperatures covered by the WNE 
and WC are similar and extend towards higher temperatures 
than the range of WNL stars. The difference in temperature 
ranges between the H-rich and the H-poor stars reflects the 
strong dependency of the opacity on the quantity of hydro- 
gen. 

Before comparing with the present evolutionary tracks, let 
us try to explain the features listed above in a general theoret- 
ical framework. We focus here on the luminosity and exclude 
the question of the effective temperature. The reason for that is 
that the effective temperature depends a lot on the mass-loss rate 
used and on the physics of the outer layers, while luminosity is 
quite tightly related to the total mass of the star and its inter- 
nal physics. In that respect, it is a more fundamental quantity to 
compare with models for the interior of stars. 

The decrease in luminosity when passing from the WNL to 
the WC stars might be interpreted in two ways. To obtain a WC 
star from a given initial mass star, the star must lose more mass 
than what is required to obtain a WNL star, implying that on 
average WC stars have lower actual masses and consequently 
also lower luminosities. Another way to interpret this feature 
(although not incompatible with the previous one) would be to 
assume that the most massive stars produce the high-luminosity 
WNL stars (which will then evolve into lower-luminosity WNE 
and WC-type stars), while the less massive stars produce the WC 
stars. In the framework of the single-star scenario, this last ex- 
planation requires very strong mass loss during the RSG phase 
of stars with initial masses of about 15 Mq. This strong mass 
loss could be caused by some physical processes originating in 
the envelope of the RSGs, while in the binary channel it could 
be caused by a RLOF process occurring in a close-binary sys- 
tem. Whatever process is invoked, it should not produce any long 
WNL or WNE phases at this low luminosity range since these 
stars are not observed. This point will be discussed in more de- 
tail in Sect. 



Comparing the above observed positions with our rotating 
stellar models we note the following features. The group of 
seven very luminous stars (log(L/L0) > 6.25) all show a mass 
fraction of hydrogen higher than 10%, except in one case where 
the mass fraction is estimated to be 0.05 (see Fig.[4]i. In addition, 
no WNE and WC are observed in the luminosity range covered 
by these stars. Is there any explanation for these two features? 
The reason why only H-rich WN stars are observed in this high- 
luminosity range supports the idea that WR stars form through 
a combination of mass-loss and radiative-zone mixing, and not 
only through mass loss as in the non-rotating models. To illus- 
trate this, we have plotted in Fig.|5]the evolution as a function of 
time of the mass fraction of hydrogen at the surface of our 60, 
85, and 120 Mq stellar models with and without rotation during 
the WNL phase. We see that only the most massive star mod- 
els with rotation (M > 60 Mq) present an H-rich surface during 
sufficiently long periods (a few 100 thousands years) for allow- 
ing this phase to be observable. Corresponding models without 
rotation, or lower initial mass models, which enter into the WR 
phase only after the end of the MS phase, do not show any long 
"H-rich" periods and thus cannot account for the most luminous 
H-rich WNL stars. The absence of WNE and WC stars in this lu- 
minosity range may be explained by the fact, already mentioned 
above, that WNE and WC stars are expected to be less massive 
and therefore less luminous than WNL stars because more mass 
has to be removed from the star to reach those stages and/or be- 
cause they are produced from stars with lower initial masses . 

Figure [4] shows the evolution of the mass fraction of hydro- 
gen at the surface as a function of the luminosity. Compared 
to non-rotating models, the rotating ones extend the regions 
covered by the WNL and WNE stars to a lower luminosity. 
However, the change remains modest so that comparisons with 
observations hardly allow one set to be favoured over the other. 
As emphasised above, the time spent in the H-rich portion of the 
diagram is probably much more decisive and favours the rotating 
models. We see that the tracks cover the region where the WNL 
stars are observed. The extension in luminosity of the WNE stars 
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Table 2. Mass ranges for the various WR- and SN-remnant-types (in Mq) deduced from our models, and comparison with previous 
works. 







O- 


ype 


WNL 


WNE 


WC 




This work 




















rot. 


15.8 


-20.0 


20.0- 


-25.3 


25.3- 


-27.0 


27.0- 


- 120.0 




no rot. 


15.0 


-25.0 


25.0- 


-31.7 


31.7- 


-40.5 


40.5- 


- 120.0 




This work (High WNE' ) 

rot. 


15.8 


-20.0 


20.0- 


-21.6 


21.6- 


-27.0 


27.0- 


- 120.0 




no rot. 


15.0 


-25.0 


25.0- 


-29.2 


29.2- 


-40.5 


40.5- 


- 120.0 




Georgv et al. (2009} 

rot. 






23.0- 


-26.0 


26.0- 


-29.0 


29.0- 


- 120.0 






SN II-P 


SN II-L/n 


SNIb 


SNIc 




i This work (low SN Ic 2 ) 

rot. 


8.0- 


•16.8 


16.8- 


-25.0 


25.0- 
39.1- 


-31.1 
120.0 


31.1 


-39.1 


> 


no rot. 


8.0- 


-19.0 


19.0- 


-32.0 


32.0- 


120.0 






05 

B/J 


This work (medium SN Ic 2 ) 
rot. 


8.0- 


- 16.8 


16.8- 


-25.0 


25.0- 


-30.1 


30.1 


-82.0 


' C 1 












82.0- 


-88.7 


88.7- 


- 120.0 


T 


no rot. 


8.0- 


-19.0 


19.0- 


-32.0 


32.0- 


-52.2 


52.2- 


- 106.4 


X 












106.4- 


- 120.0 






CO 


This work (high SN Ic 2 ) 

rot. 


8.0- 


-16.8 


16.8- 


-25.0 


25.0- 


-29.0 


29.0- 


- 120.0 




no rot. 


8.0- 


-19.0 


19.0- 


-32.0 


32.0- 


-44.8 


44.8- 


120.0 


z 

05 


r This work 


















no rot. 


8.0- 


-19.0 


19.0- 


-32.0 


32.0- 


-43.8 






BfJ 


This work (low SN Ic 2 ) 

rot. 


8.0- 


-16.8 


16.8- 


-25.0 


25.0- 


-31.1 


31.1 


-33.9 




C 

T 

03 


This work (medium SN Ic 2 ) 
rot. 


8.0- 


-16.8 


16.8- 


-25.0 


25.0- 


-30.1 


30.1 


-33.9 


This work (high SN Ic 2 ) 

rot. 


8.0- 


-16.8 


16.8- 


-25.0 


25.0- 


-29.0 


29.0 


-33.9 




Georav et al. (2009} 

rot. 




8.0- 


25.0 




25.0- 


-39.0 


39.0- 


120.0 



Notes. All masses are given in Mq. 

(1) Considering that the surface He abundance limit between WNL and WNE is 0.1 instead of 10~ 5 . 

(2) The maximum He mass ejected allowed to be still considered as a Type Ic SN is 0.4 - 0.6 - 0.8 Mq in the low - medium - high SN Ic case 
respectively (see Georav et al. 2009). 

(3) Assuming that the formation of a BH during the collapse has no influence on the SN explosion. 

(4) Assuming that the formation of a BH during the collapse prevents a bright SN explosion. 



is also well reproduced, somewhat supporting the present single- 
star models for explaining these populations. 

A serious difficulty is the low observed luminosity of some 
WC stars. Our present tracks predict a lower luminosity limit for 
the WC stars of 5.35 (in log(L/LQ)), while the lowest luminosity 
plotted in Fig. [3] is about 4.9 a ccording to the rev ised spectral 
analyses of galactic WC stars bv lSander et al.l(l2012l) . We discuss 
this point in more detail below. 



4.2. Discussion of various possible origins for the 
low-luminosity WC stars 

What might be the origin of the observed low-luminosity (and 
thus low-mass) WC stars? If the luminosity is not underesti- 
mated at present, we can imagine three different kinds of evo- 
lutionary scenarios: 



1. Stars with initial masses between 15 and 20 Mq lose much 
more mass than in the present grid of models during their 
RSG stage. We shall refer to this scenario as the RSG sce- 
nario. 



2. Stars with initial masses above about 25 Mq lose much more 
mass than in the present grid of models (the massive star 
scenario). 

3. The low-mass WC stars are produced in close-binary sys- 
tems through RLOF. This process can indeed produce lower 
final masses, as il l ustrate d, for instance, by the recent mod- 
els of lYoon et al.l yOlO). We shall refer to this scenario as 
the close-binary (CB) scenario. 

Let us now discuss the advantages and disadvantages of each of 
these three scenarios. 

1. The RSG scenario. There are some mass-loss rate de- 
terminatio ns for RSG pointin g towards very high values. For 
instance, Ivan Loon et al.l (120051) obtained that dust-enshrouded 
RSGs present mass-loss rates that are a factor of 3-50 t imes 
higher than the rate of de Jager et al. (1988). Humphrey s et al.l 
dl997ll2005l) and lSmith et alJd2009h also indicated that the RSG 
star VY CMa went t hrough important ep isodic mass ejections 
500-1000 years ago. iMoriva et al.l (1201 lb reached similar con- 
clusions, studying the luminosity curve of SNe with RSG pro- 
genitors. Their results indicate that some mechanism is proba- 
bly inducing extensive mass loss (greater th an lO^Mfriyr" 1 ) in 



massive RSGs just before their explosions. Humphreys (2008) 
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Fig. 3. Po sitions of obs e rved W N and WC stars in the HRD as 
given by Hamann et al. (2006) and Sander et al. (2012), respec- 
tively. The empty circles are WNL stars, and the full triangles 
are WNE. The WC stars are represented by pentagons, filled 
when the distance is known and empty when it is unknown. The 
present rotating tracks are superposed. 
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Fig. 4. Evolution of the mass fraction of hydrogen at the sur- 
face as a function of the luminosity. The continuous lines are 
the present rotating models, while the dotted lines represent the 
non-rotating ones. The dots are the WN stars with non-zero 
H-abundance and the shaded zone shows the range in lumi- 
nosity of the WN s tars with no H detected at the surface by 
Hamann et al. (200 6]). Th e triangles are LBV stars dGroh et all 
2009; Lamers e t al.ll200ll) . In this graph, models start at the top 
on the MS. They evolve to the right as their luminosity increases 
and downwards with mass loss peeling off the hydrogen-rich lay- 
ers. The most massive models evolve back to the left (decreasing 
luminosity) due to the strong mass loss in the WR phase. 




3 3.5 

time [Myr] 

Fig. 5. Evolution of the mass fraction of hydrogen at the surface 
as a function of time during the WNL phase. The continuous red 
lines correspond to rotating models, while the dashed black lines 
correspond to non-rotating ones. The corresponding masses are 
indicated near the curves. 



suggests that convective/magnetic activity may be the cause for 
what does appear as episodic and localise d mass-loss events . 
Also very interesting are the observations bv lDavieset al. ( 2008) 
of the so-called RSG clusters in the direction of Scutum. Maser 
emission indicative of a high-density medium above the photo- 
sphere is observed around the most luminous RSGs. This can be 
caused by the high mass-loss rates experienced by these stars. 
The presence of a luminous yellow supergiant in one of these 
clusters is also consistent with the idea that this star evolved 
away from the RSG stage (Davies et al. 2008). 

iMauron & Josse lin (201 lL on the other hand, still recom- 
mend using the Ide Jager et a l. (1988) rate for Galactic RSGs, 
indicating that this prescription agrees to within a factor of 4 
with most mass-loss rate estimates based on the infrared 60 /im 
excess. This result is however compatible with the existence of 
short phases during which the mass-loss rates are much stronger. 
Indeed, mass-loss rates during the RSG phase may present some 
outbursting characteristics, somewhat similar (at lower luminos- 
ity and effective temperature) to what happens at higher lumi- 
nosity and e ffective temperatur e for the luminous blue variable 
(LBV) stars dSmith et aLll201 lb . The physical reasons for these 
high mass-loss rates may be re lated to the pulsational proper- 
ties of RSGs dYoon & Cantielldl20 10) and/or to the appearance 
of super -Eddington luminosities in the outer layers of the star 
(iPanerJ) . 

From the theoretical point of view, it is well known that 
strong mass loss during the RSG phase favours a bluewards evo- 
lution dSalasnichetal.l ll999; Van beveren et al.| [T998. 2007) and 
thus helps in fulfilling one of the two minimal conditions for hav- 
ing a WR star, namely to have an effective temperature higher 
than about log(r e ff/K) = 4.0. However, it is unclear whe ther this 
strong mass loss may lead to the formation of WC stars. iGeorgvl 
(2012) has shown that increased mass-loss rates during the RSG 
phase (up to a factor of 10) are insufficient to lead to a WR star 
at the end of the evolution. To be valid, the RSG scenario should 
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therefore involve a very strong increase of the mass-loss rate of 
more than 10 times during that phase with respect to the standard 
one. 

One can already note that it will be difficult to separate the 
physical reason for the envelope loss, determining whether it is 
caused by a physical process occurring in the envelope of the star 
(as would be the case in the single-star channel) or to a RLOF. A 
constraint on this scenario may come from the fact that there are 
very few (if any) si ngle-age clusters that si multaneously show 
RSGs and WR stars (Humphreys et al. 1985), except in the clus- 
te rs at the centre of our G alaxy (lFigerl i2007) and in Westerlund 
1 dNegueruela et al.l2010. although in the latter case we may see 
two clusters aligne d instead of only on e, as in the case of the 
Dansk clusters; see iDavies et al.ll2.Q12b . If some WR stars are 
produced from a range of initial masses that are also producing 
RSG, then to be consistent with this observation, either the RSG 
duration or the WR duration should be very short. 

2. The massi v e star scenario. In the high mass-loss rate grid 
of iMevnet et alj (11994). a very low luminosity is reached for 
WC stars due to heavy mass loss. The lowest luminosity for 
the WC stars was at 4.5 and originated fro m the evolution of a 
12OM0! The high mass-loss rates used in IMevnet et a l. (1994) 
are no longer supported by the more recent mass-loss determi- 
nations for O-type and WR stars that account for th e effects of 
clumping dVink et aUbOOlt iNugis & Lamersll2000l) . However, 
very massive stars could lose very large amounts of mass in 
very short periods durin g which a strong mass outburst occurs 
dSmith & O wocki 2006). Therefore, it is possible that the mass- 
loss determination, necessarily based on more frequent "normal- 
moderate" mass-loss rate stages, actually underestimates the true 
time-averaged mass-loss rate. 

3. The CB scenario. In close binaries, the primary can lose 
its H-rich envelope during RLOF phases. The secondary can also 
undergo such a loss. The two stars may also merge or enter into 
a common-envelope phase leading to heavy mass loss episodes. 
Clearly, these evolutionary scenarios can lead to the produc- 
tion of WR stars and therefore have an impact on their pop- 
ulation s. For instance, the close-binary scenario of Yoo n et alj 
(2010) makes it possible to produce final masses in the range be- 
tween 1 and 1 Mq, even when starting from hig h initial masses 



(as for instance a 60 M ^, see Yoon et al.l 12010). However, ac 



cording to the review by ICrowfherl (120071 see Fig. 4), the least 



massive WC star whose mass has been determined from binary 
orbit has a mass of about 9 Mq (the most massive has about 
16 Mq). The recent analysis of Galactic WC stars performed by 
Sander et al. (2012) gives masses between 8 and 31 Mq, cov- 
ering quite well the masses obtained by our single-star models 
(10-26 Mq). Therefore, the least massive He-core produced by 
the close-binary scenario cannot be invoked to explain the low- 
luminosity WC stars. 

Close-binary evolution can produce WR stars from lower 
initial ma sses than the single-star evolutionary channel. For in- 
stance, Eldrid ge et al.l (12008) found that the minimum initial 
masses for stars to become WC stars is lowered from about 
27 Mq in the single-star scenario to 15 Mq in the close-binary 
scenario. In the mass range between 15 and 25 Mq, the WC 
stage occurs only at the very end of the evolution and has a du- 
ration of the order of 10 4 years. This supports the view that low- 
luminosity WC stars could indeed be the result of close-binary 
evolution. It remains to be seen whether this scenario is able to 
explain the observed number of low-luminosity WC stars as well 
as their masses. 

Accordingly at the moment, the three scenarios described 
above may contribute in explaining the low-luminosity WC 



stars. To investigate the relative weight of these different scenar- 
ios, progress must be made in several directions. We need to con- 
strain the occurrence of short and strong mass-ejection episodes 
during the RSG phase. Could these events allow stars with initial 
masses around 15 Mq to evolve into a WR phase? 

Is there any evidence that the low-luminosity WC stars are 
produced in close-binary systems? Positive evidence could in- 
clude the presence of a companion, or the determination that 
these stars are runaway stars, kicked off when the primary ex- 
ploded in an SN event. Note however, that runaway stars may 
also be produced by dynamical interactions in dense clusters and 
thus result from processes other than binarity. 

Is there any evidence for the presence of low-luminosity 
WC stars in very young associations, typically with a mass at 
the turnoff above about 25 MqI These low-luminosity WC stars 
would originate from an initial mass more massive than 25 Mq 
and this would support the second scenario (the massive star sce- 
nario). 

4.3. Surface chemical compositions 

In Fig. [6] the evolution of the surface abundances is shown 
for the n on-rotating and rot a ting 6 Mq mod els. As already 
noted bv IMevnet & Maeder (120031) (see also iMaederi 11987b 
iFliegner & Langerl 19951) . one of the main effects of rotation is to 
smooth the internal chemical gradients and to favour a more pro- 
gressive arrival of internal nuclear products at the surface. This 
is well apparent comparing the curves shown in the left and right 
panels of Fig. [6] 

One observable consequence is that rotating models predict 
more extended phases during which both H- and He-burning 
products are seen at the surface ( Langerlll991b . Note that in the 
left panel of Fig. [6] the WNC phase is so narrow that it only ap- 
pears as a thicker tick in the upper part of the top panel, while in 
the right panel of Fig. [6] a well-extended phase is present. This 
phase overlaps with the WNE phase, which is not shown in the 
figure for clarity. 

Except for this effect, rotation leaves no other easily observ- 
able imprints on the way the surface composition evolves. This 
is quite expected because CNO equilibrium values are obtained 
during the WN phase that can be deduced from the nuclear prop- 
erties of the chemical species and are not strongly affected by the 
details of the considered stellar model. 

This is also true for the abundance of 22 Ne obtained during 
the WC phase. The abundance of 22 Ne at this stage (WC) is actu- 
ally an indication of the initial CNO content of the star. Indeed, 
22 Ne comes mainly from the destruction of 14 N at the beginning 
of the helium-burning phase. This 14 N is the result of the trans- 
formation of carbon and oxygen into nitrogen operated by the 
CNO cycle during the core H-burning phase. 

It is interesting to note that comparing the observed Ne/He 
ratios at the surface of WC stars with models computed with 
Z = 0.020 shows that models over-predict the Ne abun- 
da nce, while models sta rting with the solar abundances given 
bv lAsplundet al.l Q2P05) yield a much better fit, as can be seen 
in Fig. [7] This confirms that massive stars in the s olar neighbour- 
hood have initial metallicities that agree with the lAsplund et al.l 
(2005) solar abundances. 

Let us note that this overabundance of 22 Ne at the surface of 
WC stars is not only an important confirmation of the nuclear 
reaction chains occurring during He-burning, but is also related 
to the questions of t he origin of the m aterial accelerated into 
galactic cosmic rays dBinns et al.l l2005) and to that of the weak 
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Fig. 6. Top left panel: Evolution of the surface abundances in mass fraction as a function of the actual mass of the star for the non- 
rotating 60 Mq model. Different evolutionary phases are indicated in the upper part of the figure. A very narrow region, just before 
the WC phase, corresponds to the WNC phase (not labelled in the figure for clarity). Bottom left panel: Evolution of abundance 
ratios (in number) as a function of the actual mass of the star. Right panels: Same as the left panels for the rotating model. 



s-process in massive stars, since 22 Ne is the main neutron source 
in these stars dPetersll 19681) . 

Rotation increases the total quantity of ejected 26 A1, as can 
be seen comparing the non-rotating and the rotating models. 
Note however, that the quantity of 26 A1 ejected by the winds of 
massive stars i s lower in the present study than that found by 
iPalacios et alj d2005bl) for a given initial mass. This is a natural 
consequence of the lower "solar" metallicity. Since the new so- 
lar abundances initially have fewer metals, this means that there 
is less 25 Mg, the seed element for the synthesis of 26 A1. The im- 
pact on the global production of 26 Al in the Milky Way therefore 
needs to be revised. This will be done in a future paper. For now, 
we can say that the revised solar abundances diminish the role 
played by the WR stellar winds in the enrichment of the inter- 
stellar medium in 26 A1. 

As was proposed by iMevnet & Arnouldl ([2000), massive 
stars can also enrich the interstellar medium in fluorine through 
their winds. We con firm this with the present models (see also 
IPalacios et alJl2005ah . Rotation does not appear to have a great 
impact on the wind ejection of this element. 

4.4. WR populations 

It is important to realise that to understand the origin of the WR 
populations, it is necessary to account for not only the relative 
number of WR of different subtypes with respect to the total 
number of WR, but also for the number of WR to O-type stars. 
For instance, if a good fit is obtained for the ratio WC/WR but 
the models only account for 10% of the observed WR/O ratio, 



then this implies that the models predict a good WC/WR ratio 
in about 10% of the WR stars with no guarantee that the fit will 
be good for the remaining 90% of the WR stars. Therefore, it is 
very important to look at as many number ratios as possible to 
assess the quality of a model. 

For the observed number ratio of WR to O-type stars, we 
focus on the spherical region of 2.5 kpc radius centred on the 
Sun for which extensive data are found in the literature. In 
the se venth catalogue of Galactic WR stars of van der Hucht 
(2001), 64 WR stars were observed in that regi o n. Ten years 
later, new detections by Pasquali et al. (2002), Negueruela 
(2003), Hopewell et al. (2005)' Hadfield & Crowther"(2007), 
Sharaetal. (2009) and Mauerhan et al. (2011) increased that 
number by about 20% to 77. In that same zone the number of 
O-type stars is given to be equal to 48 1 in the iGarmanv et al.l 
(1982) catalogue. However, there is evidence that this number 
significantly underestimates the number of early-type stars in 
that region. For instance, more recent observations analysed by 
iReedl (1200 1) indicate that the number of OB stars more massive 
than 20 Mq within 2. 5 kpc of the Sun might be 50% greater than 
the value obtained by Garma nv et al.l (11982). In Table[3]we use a 
number of O-type star within 2.5 kpc from the Sun increased b 
33% with respect to the number given by G armanv et 



:reased by 
ID ([1982) 



to account for some incompleteness. This gives a WR/O number 
ratio equal to 0. 12 (=77/636). An error of about 25% around this 
value is probably reasonable. 

For the WR subtype populations, we consider the 3 kpc ra- 
dius region centred on the Sun. Using the data collected in the 
above references, plus a few recent detections in the region be- 
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Fig. 7. Variations of the abundance ratios Ne/He vs. C/He at the 
surface of WC stars (in number). Models are represented by the 
lines (the Z=0.014 are the rotating models of the present paper, 
the Z = 0.020 is the rotating model of Mevnet & Maeder 2003). 
The points are the observed v alues bv lDessart et al.ld2000l) (filled 
circles) and lCrowtherl (120061) . 



tween 2.5 and 3 kpc (Mauerhan et al. 2009; Roman-Lopes 201 1; 
lAnderson et al.ll201 II) . there are 45 WN and 55 WC stars as dif- 
ferent subtypes of WR stars. The 45 WN stars are distributed 
between 21 WNL, 23 WNE and 1 WNC stars. It must be em- 
phasised here that the respective portion of WNL and WNE stars 
is quite uncertain. Where information on the mass fraction of H 
was available, we used the presence of H to classify a star as a 
WNL star and the non-detection of hydrogen for classifying it as 
a WNE star. For this we us ed the sample of Galactic WN stars 
analysed by Hama nri et al.l d 20061) . which has 20 stars in com- 
mon with our sample of nearby WR stars. We classified in the 
remaining 25 stars all stars with an "h", indicative of the pres- 
ence of hydrogen, as WNL, as well as all stars with spectral type 
equal or later than WN7; the remaining stars are considered to 
be WNE stars. Again, assuming an error of about 25% we obtain 
the range of values given in Table [3] The WNC star is special in 
the sense that there is only one star of this type in the solar neigh- 
bourhood, which represents 1% of WR stars in that same region. 
Among th e WR stars presently d etected in the Galaxy (476 ac- 
cording to Mauerh an et al. 20111 b ut 548 if we add the 72 new 
WR detections by Shara et al .1120 lib . 9 are at present identified 
as WNC stars, which corresponds to a proportion of about 2%. 
Therefore we indicate a range between 1 and 2% for this ratio in 
Tabled 

To compute the ratio of two types of stellar populations (A 
and B), we have to know the initial mass range leading to each 
population [m a ...M a .„1 and |m b . ...Afl„l, as well as the 

" " L mm maxj [ mm maxj' 

time spent in the corresponding phase r A and r B . When a con- 
stant star formation rate is assumed, the ratio of the two popula- 
tions is given by 



N(A) _ Im£ <P{M)rHM)dM 
N < B > = £f- >(M)rB(M)dM' 



where </>(M) is the initial mass funct ion. Here w e consider a 
Salpeter initial mass function (IMF) (Salpet eril955l) . The results 
are presented in Table [3] 

The predicted value for the rotating models corresponds to 
only one initial Vi n i/v cr i t value (0.4). A detailed computation 
should account for a distribution of the initial velocities. Since 
the value of 0.4 has been chosen to reproduce the average ob- 
serve d surface velocities during the MS for B stars (IHuang et al.l 
2010), we expect that the values of the ratios obtained for this 
particular value are close to the value that would be obtained by 
a properly weighted averaging over an initial velocity distribu- 
tion (assuming that the characteristic initial ratio Vj n i/v cr i t is the 
same for more massive stars). 

The value of the WR/O ratio given by the non-rotating 
single-star model is 0.02, which is quite low compared to the 
observed ratio of around 0.12. The low ratio results from the 
low mass-loss rates used here that account for the clumping ef- 
fect. If these models were the correct ones, then it would mean 
that close-binary evolution is responsible for more than 80% of 
the WR stars. However, since rotating models better fit many 
observed features of massive stars (such as, for example, the 
change of the surface abundances), their predictions are to be 
preferred to those obtained from the non-rotating models. The 
WR/O value obtained by the rotating models, 0.07, is more than 
three times the value obtained by the non-rotating models. It is 
below the observed values, which leaves some room for ~ 40% 
of binaries from the close-binary scenario to contribute to the 
number of WR stars at solar metallicity. 

Since the WNL phase is increased in duration by rotation, the 
WNL/WR ratio is increased when the rotating models are used. 
The WNE/WR and WC/WR ratios, in contrast, decrease. The 
few WNC stars can be very well explained by internal mixing 
processes. 

To confront the present theoretical predictions with the 
observed numbers , we followed the same formalism as in 
iMaeder & Mevnetl (11994 . WR S is the number of WR stars that 
do not require any mass transfer episode for entering into their 
WR phase, and WR c b are those WR stars that owe their WR na- 
ture to a previous RLOF mass transfer. The total ratio WR/O can 
be written as the sum WR s /0 + WR cb /0. <p is the ratio WR cb /0. 
Below we consider that the WNC stars are produced by internal 
mixing only, not by a binary mass-transfer event. In principle, a 
WNC star might result from a WN star that accretes matter from 
a WC star, but this does not appear to be realistic for at least two 
reasons: first, because of the compactness of these types of stars, 
some very peculiar initial conditions for the mass ratio and the 
period of the orbit are likely necessary; second, both types of 
stars suffer strong stellar winds, and processes such as wind col- 
lisions are more likely to occur than accretion. The WNC/WR 
ratio will therefore depend on <p only through the dependence of 
the total number of WR. 

The first three upper panels of Fig. [8] show how the WR/O, 
WNC/WR and WR c b/ WR0 compare to observations for different 
values of ip, using the results of the present rotating models for 
the fractions obtained when ip - 0. We see that the range of 
values for ip compatible with the observations is between 0.03 
and 0.075, which corresponds to situations when 31 to 53% of 
the WR stars owe their WR nature to a RLOF episode. 

To further pursue this line of reasoning and use comparisons 
with the observed WNL, WNE and WC star number ratios, we 
need to make some assumptions about the distributions in the 



(1) 



3 The fraction of WR due to close-binary evolution with respect to 
the total number of WR stars is given by <p/(WR s /0 + tp). 
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0.3 



N(WNE)/N(WR)- 

H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




Fig. 8. Variations of various number ratios as a function of 
ip - WR c b/0, the fraction of WR stars with respect to O-type 
stars that owe their nature as WR stars to a RLOF episode in a 
close-binary system. The continuous lines in the 4 lower panels 
are the models for y?wNL, Vwne, Vwnc and ipwc equal to 0, 0, 
Q3(p, and OJtp respectively (see text). The horizontal strips cor- 
respond to observed ratios as reported in Table [3] The green part 
corresponds to the values of <p allowed by the observation, and 
the red one the excluded values. 



WR subtypes resulting from the close-binary channel. We de- 
fine i^wnl as WNL c b/WR c b. Analog definitions are considered 
for the other WR subtypes. With these definitions, the ratio of 
WR(st)/WR where "st" designates a given subtype can be writ- 
ten 

WR s (st) 

WR(st) _ ^mr 



WR 



1 



wr s /o 



WRJO 



The values of the different components of <p should of course be 
the outcome of close-binary e v oluti on calculations, such as those 
performed by Eldridge et al. (2008). Here, however, let us try 
a simpler approach guided by the comparisons of the observed 
ratios with the predictions of the rotating models. 

For the reasons laid out above, we assume that <^>wnc = 0. 
The WNL/WR ratio given by the present single-star rotating 
models is well above the observed range. This ratio therefore 
needs to be reduced when the close-binary scenarios are ac- 
counted for. To go in that direction we make the extreme as- 
sumption that i£>wnl = 0. Even so, we see that a very high frac- 
tion of binaries would still be needed, which would be clearly 
incompatible with the observed WR/O. This may indicate that 




N(WNE)/N(WR) 

H 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 




Fig. 9. Same as in Fig. [8] but considering that WNE stars are 
WN stars with less than 10% of hydrogen at the surface (mass 
fraction). 



the WNL phase is somewhat overestimated in the present mod- 
elfl 

Considering that </?wne = 0.3^ and i^wc = 0.7y, we find a 
consistent solution for if - 0.075 for all ratios shown in Fig. [8] 
with the exception of the WNL/WR ratio, which is above the 
observed range. At this point we conclude that mainly because 
of the high WNL/WR ratio obtained from the present models, 
we do not succeed in finding a completely consistent solution. 
To improve the situation, two lines of investigations can be fol- 
lowed, one observational and one theoretical. 

From the observational point of view, it is difficult to match 
theoretical criteria and observational ones. For example, how 
much H can be hidden in an observed WNE? Obtaining new reli- 
able measurements of the H abundance in all WN stars of the so- 
lar neighbourhood is crucial to find a better criterion for linking 
a given structure obtained in stellar models to an observed WR 
subtype. Some stars have changed classifications when more re- 
fined observations have been performed, indicating that the dis- 
tinction between WNL and WNE is not always clear-cut. To il- 
lustrate how the WNL/WR and WNE/WR fractions depend on 
the definition of WNE, we indicate in Table[3]the ratios obtained 
assuming that WNE stars are all WN stars with less than 10% 
of hydrogen at the surface ("high- WNE case", see Tables |2]and 
[3). We also show in Fig. [9] how the theoretical ratios compare 



4 We recall that here we only consider one value of v ini /v crit ; the 
value obtained by a proper averaging over an initial velocity distribu- 
tion would likely be lower than that obtained from the unique value of 
Vim/Vcrit = 0.4 
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with the observations. When adopting this new definition for the 
WNE stars, the WNE/WR ratio increases from 0.022 to 0.147, 
nearly a factor of 7. A consistent solution is now obtained for all 
ratios and for if ~ 0.07, which corresponds to a share of about 
50-50 between the single and binary channel for producing the 
observed WR populations at Z = 0.014. 

From the theoretical approach, one would explore a different 
mass-loss rate history during the WNL and WNE phases. For 
instance, are the WN mass-loss rates the same when the star is a 
core H-burning star or a core He-burning one? Are the mass-loss 
rates during the WNE phase overestimated at present? Stellar 
models in the future may explore the consequences of different 
prescriptions and offer some guidelines for improving the situa- 
tions. 



5. Progenitors of type Ibc supernovae 

The questions that we shall address in this section are the fol- 
lowing ones: 

- What are the predictions of the present models for the fre- 
quency of type lb and Ic SNe? 

- How do these predictions compare to the observed ratios? 

- Can we find a consistent picture for understanding both the 
WR populations as observed in the solar neighbourhood and 
the frequency of type Ibc SNe, including contributions of 
close-binary evolution? 

- Which fraction of the present models would produce a 
fast rotating BH, accompanied by a type Ic SNe? The last 
question is triggered by the possibility that stars showing 
these condition s may give birth to a collapsar, proposed by 
IWooslevi (1 19931) as progenitors of LGRBs. 

5.1. Predictions of the models 

To determine the type of the SN that occurs at the end of the 
evolution of o ur massive star mo dels, we used the same pro- 
cedure as in IGeorgy et a l. (2009). We recall here a few main 
points. First, we obtain the baryonic mass of the remnant us- 
ing the same relation betw een the CO-cor e mass and the rem- 
nant baryonic mass as in iMaederl (U992). This mass is then 
used to compute th e gravitationa l mas s of the remnant, with 
the relation given in Hirschi et al. (2005). The results are given 
in Table [4] Using these data, we estimate the maximal mass 
on the zero-age main se quence (ZAM S ) pro ducing an NS dur- 
ing the SN event. As in Ge orgy et al.l J2009I) . we consider the 
most massive NS to be 2.7 Mq. This choice is somewhat sup- 
ported by the recen t disco very of a very massive 2.4 Mq NS by 
Ivan Kerkwiik et alJ (l201ll) . We find that M^ BH = 33.9 Mq (ro- 
tation), and M,^ BH = 43.8 M© (no rotation). The difference is 
explained by the rotational mixing, which increases the size of 
the core. These limits are only poorly constrained, however, due 
to uncertainties during the explosion process, which are not ac- 
counted for in our simple estimate (for example, the amount of 
matter that falls back, etc., see Frver 2006). 

Once the baryonic mass of the remnant is known, we deduce 
the composition of the ejecta, assuming that the entire mass be- 
tween the surface and the edge of the remnant is ejected during 
the SN event. The results are shown in Table|4] In the same table, 
we also indicate the SN type, using the chemical composition 
of the ejecta and following the same criteria as in IGeorgy et al.l 
J2009I) : 

- If there is some hydrogen in the ejecta, the SN is a type II. 



- If there is no hydrogen and a helium mass higher than a given 
value (M^), it is a type lb. 

- If there is no hydrogen and a helium mass lower than MJft, 
it is a type Ic. 

Accounting for the uncertainty of MS, we give three values for 
the mass limit between types lb and Ic SNe, with Mj^ — 0.4 Mq 
(low SN Ic, as this criterion favours type lb SNe), M^ = 0.6 Mq 
(medium case), and M^ = 0.8 Mq (high SN Ic, as it favours 
type Ic SNe). 

Recently, iDessart et al.l (|2011|) found that some amount of 
He can be unobservable in the spectrum of an SN when it is lo- 
cated primarily in the most external layers (up to 50% of He in 
mass fraction in the first 1 Mq immediately below the surface). 
In this case, our non-rotating 120 Mq and rotating 85 Mq mod- 
els would be classified as a type Ic SN instead of type lb. This 
criterion thus produces the same mass limit between type lb and 
type Ic SNe as in the "high SN Ic" case (see Tabled. 

Contrarily to Georgv et al. (2009), where no distinction was 
considered between the various subtyp es of type II SN e, we add 
in this work the same criterion as in lHeger et al.l (12003) to distin- 
guish between the type II-P SNe, and the type II-L or Il-b SNe 
(owing to the weak statistics of these events, there is no simple 
way to distinguish between type II-L and II -b on the basis of the 
ejecta composition only. We therefore consider them as a unique 
sample): 

- A type II SN is considered to be a type II-P if the ejecta 
contains more than 2 Mq of H. 

- In the other case, it is considered to be a type II-L/b. 

We see from Table |4] that the upper mass limit for type II-P SN 
would be between 15 and 20 Mq for both non-rotating and rotat- 
ing models. More precise mass limits can be obtained by inter- 
polation. The mass limits become 19.0 Mq for the non-rotating 
models and 16.8 Mp, for the ro tating ones. Both values agree 
with the results bv lSmartt etall (l2009h . 

We computed the var i ous S N type ratios with the same 
method as in Georgy et al. (2009), with the same IMF. We con- 
sidered two cases: 1) a SN event is visible even when a BH is 
formed; 2) no visible SN occurs when a BH is formed. The re- 
sults are shown in Table [3] Th ey are quite similar to those ob - 
tained with different models bv lMevnet & Maedej(l2003ll2005l) . 
IGeorgy et al. (2009). Therefore, these results show some robust- 
ness against many changes in the physical ingredients of the 
models. 

We see that rotation increases the SN Ibc/SN II ratio by about 
60% (48%) in case 1 (2). The share of SN type Ibc between the 
lb and the Ic types depends on the value chosen for M^. A 
low value implies more restrictive possibilities to obtain a SN 
Ic and consequently a lower SN Ic/SN II ratio and a higher SN 
Ib/SN II ratio. The trend of having more Ic SNe when rotation is 
accounted for remains true for all M, He considered here: rotation 

hm 

decreases the minimum initial mass of stars finishing their life 
as WC (see Table 0. The effect of rotation on the SN Ib/SN II 
ratio is less clear: the ratio increases or decreases when rotation 
is accounted for depending on Mjfe, 



5.2. Comparisons with observations 

At solar metallicity, in case 1, the present models give ratios 
well within the range of the observed values determined by 
iBoissier & Prantzosl (120091) . Indeed, one sees that the present 
single-star models may explain the greatest part of the observed 
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Table 3. WR- and SN-type ratios deduced from our models, and comparison with previous works and observations. 







WR/O-stars 


WNL/WR 


WNE/WR 


WN/WR 


WNC/WR 


WC/WR 


WC/WN 




This work 


















rot. 


0.066 


0.687 


0.022 


0.709 


0.028 


0.291 


0.409 




no rot. 


0.015 


0.253 


0.168 


0.421 


0.003 


0.579 


1.376 




This work (High WNE) 1 


















rot. 


0.065 


0.562 


0.147 


0.709 


0.028 


0.291 


0.410 




no rot. 


0.015 


0.185 


0.238 


0.423 


0.003 


0.577 


1.364 




Mevnet & Maeder (2003) 


















rot. 


0.07 


0.66 


0.05 


0.71 


0.04 


0.25 


0.35 




no rot. 


0.02 


0.35 


0.16 


0.51 


0.00 


0.49 


0.97 




Observations 


















Solar Neighbourhood 


0.12±0.03 


0.21±0.05 


0.23±0.06 


0.44±0.11 


0.015±0.05 


0.55±0.14 


1.25±0.30 






SN Ibc / SN II 


SN lb / SN II 


SN Ic / SN II 


SN Ic / SN lb 










' This work (low SN Ic 2 ) 












m 


rot. 


0.241 


0.186 


0.054 


0.290 








z 

as 


no rot. 


0.151 


0.151 


- 


- 








Jc 


This work (medium SN Ic 2 ) 
















■~ ' 


rot. 


0.241 


0.066 


0.174 


2.635 








t 

X 


no rot. 


0.151 


0.093 


0.058 


0.622 








This work (high SN Ic 2 ) 
















n 


rot. 


0.241 


0.050 


0.191 


3.845 










^ no rot. 


0.151 


0.066 


0.085 


1.279 








Z 


This work 
















no rot. 


0.062 


0.062 


- 


- 








Jc 


This work (low SN Ic 2 ) 
















■" 


rot. 


0.092 


0.070 


0.022 


0.320 








c 


This work (medium SN Ic 2 ) 
















t 

x 


rot. 


0.092 


0.061 


0.031 


0.521 








This work (high SN Ic 2 ) 
















CD 


rot. 


0.092 


0.050 


0.042 


0.857 










Georev et al. (20091 


















rot. 


0.241 


0.060 


0.181 


3.032 










Mevnet & Maeder (20031 


















rot. 


0.32 
















no rot. 


0.12 
















Observations 


















Boissier & Prantzos (20091 


0.14-0.26 






0.39-1.13 










Prieto et al. (20081 (Z ~ 1.3 Zm) 


0.572 ± 0.286 
















Prantzos & Boissier (2003) 


0.14-0.33 















Notes. (1) Considering that the surface He abundance limit between WNL and WNE is 0.1 instead of 1(T 5 . 

(2) The maximu m He mass ejected allowed to be still considered as a type Ic SN is 0.4 - 0.6 - 0.8 A/q in the low - medium - high SN Ic case, 
respectively (see lGeorev et al.ll2009h . 

(3) Assuming that the formation of a BH during the collapse has no influence on the SN explosion. 

(4) Assuming that the formation of a BH during the collapse prevents a bright SN explosion. 

Remark: note that even though it seems that the non-rotating models better reproduce the various ratios of the WR subtypes, they are worse at 
reproducing the WR/O ratio. It is important to fit all these ratios simultaneously. 



number of type Ibc (normalised to the number of type II SNe). 
For the SN Ic/SN lb ratio, as already noted above, the value 
greatly depends on the value of M|ft, We can only say that to 
achieve a reasonable fit, rotating models would favour a low- 
medium value while the non-rotating ones favour a medium-high 
value. 

In our case 2, about half of the type Ibc SNe occurring at so- 
lar metallicity could originate from the singl e-star channel (sim- 
ilar to the conclusion we reached in iGeorgv et al.l 12009). The 
majority of the type Ibc SNe produced by the single-star chan- 
nel would be type lb SNe. If the hypothesis "formation of BH 
implies no visible SN" is correct, a significant fraction of the Ibc 
SNe must be produced in scenarios different from those com- 
puted in the present grid of models. They must originate from 
stars producing sufficiently low-mass cores to avoid BH forma- 
tion. Note also that if the "massive star" scenario is correct (see 
Sect. 14.2b . this would lead to less massive cores, and thus favour 
the production of NS instead of BH. 

The idea of very low-luminosity pro genitors for ty pe Ibc SNe 
has recently received some support by ISmartn (120091) . This re- 



view highlights that no progenitors have been detected at the 
positions of ten type Ibc SNe for which images have been ob- 
tained prior to the SN event. As indicated by this study, such 
a non-detection cannot rule out a massive WR star progenitor. 
However, under the hypothesis that the progenitor population of 
all Ibc SNe are massive WR stars, the probability that they did 
not detect any of the 10 progenitors by chance would then be 
only 11%, assuming that the lb progenitors are WN stars and I c 
progenitors are WC/WO stars (see more detail in ISmarttl 12009). 
Of course, uncertainties about the distance and the estimated de- 
gree of absorption (owing, for instance, to an enhanced mass 
loss episode just before explosion) may play a role when reject- 
ing typical WR stars as good candidates for the progenitors of 
type Ibc SNe. 

Let us suppose, for discussion purposes, that at least some 
of the progenitors of type Ibc SNe could be low-mass, low- 
luminosity, naked CO-core stars. We may first wonder to which 
extent such a population is observed. Could it correspond to the 
low-luminosity WC stars discussed above? The answer is likely 
no, since many of the low-luminosity WC stars analysed by 
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Table 4. Properties of the models at the end of their evolution. 



Mini 


Vini /Vcrit 


M He ' 


M co 


M rem (baryon./grav.) 


Tot. Ej. 


H 


He 


Prog, type 


SN type 2 


Remnant 


Spills 


M e 




M 


M 


M 


Mq 


Mq 


M 








s 


120 





(30.91) 


30.13 


9.13/4.51 


21.72 


0.00 


0.67 


WC 


lb 


BH 


- 




0.4 


(19.04) 


18.46 


5.69 / 3.65 


13.29 


0.00 


0.52 


WC 


Ic 


BH 


- 


85 





(18.65) 


17.98 


5.54/3.60 


13.09 


0.00 


0.47 


WC 


Ic 


BH 


- 




0.4 


(26.39) 


25.72 


7.88/4.26 


18.47 


0.00 


0.61 


WC 


lb 


BH 


- 


60 





(12.50) 


12.23 


3.91/2.91 


8.57 


0.00 


0.42 


WC 


Ic 


BH 


- 




0.4 


(17.98) 


17.52 


5.40/3.55 


12.58 


0.00 


0.50 


WC 


Ic+coll.? 


BH 


- 


40 





(12.82) 


10.08 


3.40 / 2.64 


9.41 


0.00 


0.95 


WNE 


lb 


NS 


- 




0.4 


(12.33) 


12.21 


3.90/2.90 


8.42 


0.00 


0.41 


WC 


Ic+coll.? 


BH 


- 


32 





(10.92) 


8.49 


3.02/2.42 


7.89 


0.00 


2.29 


WNE 


lb 


NS 


- 




0.4 


(10.13) 


10.01 


3.38/2.63 


6.73 


0.00 


0.28 


WC 


Ic 


NS 


1.3- 10~ 4 


25 





8.12 


5.95 


2.41 / 2.03 


5.85 


0.03 


2.20 


WNL 


II-L/b 


NS 


- 




0.4 


9.69 


7.09 


2.69 / 2.22 


6.99 


0.00 


1.59 


WC 


lb 


NS 


7.9 ■ 10~ 5 


20 





6.21 


4.00 


1.91/1.68 


6.66 


1.15 


3.31 


RSG 


II-L/b 


NS 


- 




0.4 


7.17 


4.73 


2.10/1.82 


5.06 


0.02 


1.61 


WNL 


II-L/b 


NS 


9.7 • 10~ 5 


15 





4.25 


2.41 


1.50/1.37 


11.72 


5.81 


4.82 


RSG 


II-P 


NS 


- 




0.4 


5.11 


3.19 


1.71/1.53 


9.36 


3.31 


3.75 


RSG 


II-P 


NS 


9.2 • 10~ 5 


12 





2.99 


1.75 


1.33/1.23 


9.97 


5.64 


3.75 


RSG 


II-P 


NS 


- 




0.4 


3.90 


2.34 


1.48/1.35 


8.73 


4.02 


3.31 


RSG 


II-P 


NS 


9.5 • 10~ 5 


9 





1.21 


1.20 


1.12/1.05 


7.64 


4.58 


2.87 


RSG 


II-P 


NS 


- 




0.4 


3.08 


1.64 


1.30/1.20 


7.22 


3.53 


3.03 


RSG 


II-P 


NS 


1.2 • 10~ 4 



Notes. The columns are: the initial mass (column 1), the initial velocity (column 2), the final masses of the He (column 3) and CO (column 4) 
cores, the baryonic and gravitational mass of the remnant (column 5), the total mass of the ejecta (column 6), the ejected mass of H (column 7) 
and He (column 8), the progenitor type (column 9), the SN type (column 10), the remnant type (column 11) and the pulsation period of the NS 
(column 12). 

(1) The parenthesis indicates that the value corresponds to the whole stellar mass. 

(2) The SN types given here are determined using the "medium SN Ic" criterion (MS = 0.6 Mq). The models that might lead to a LGRB through 
a collapsar event are indicated by the label "coll.?" (see text). 



Sand er etal .1 (12012) have visual magnitudes well above the lower 
vi sual magnitude , which would imply detection of the progenitor 
bylSmartt (2009). This means that the progenitors of these type 
Ibc could be completely different from the WR stars we observe, 
even at the lowest luminosity. 

Why do we not observe such progenitors? Where are they? 
They may have escaped detection in the case of the ten Ibc SNe 
because they are too distant in these particular cases, but how 
can we explain the fact that this population has escaped detec- 
tions independent of any link with particular SN events? He- or 
CO- rich stars of a few solar masses are luminous enough to be 
detected as individual stars. If they are not observed, an expla- 
nation has to be provided. Are they hidden in the light of their 
companion in a close-binary system? Do they escape detection 
because they may be enshrouded in heavy circumstellar mate- 
rial coming from their companion, or from matter lost by the 
system? Is their lifetime as a naked He- or CO-core so short 
that they are very seldom caught in that stage of their evolution? 
These questions remain quite open at the moment. 

If these stars exist, they would likely be produced in close- 
binary systems. As mentioned above, the RLOF process makes 
it possible to produce very low final masses between 1 and 
7 Mq, even when starting from high initial masses (for instance 
a 60 Mq, see Yoon et al. 2010). It also makes it possible to pro- 
duce naked He -cores from low er initial-mass stars. Close-binary 
scenarios of I Yoon et al.l (J2010) actually predict a bimodal distri- 
bution for the progenitors of type Ic: the high-mass progenitors 
(initial masses greater than 35 - 45 Mq) and the low-mass pro- 
genitors (between 12.5 - 13.5 Mq), However, this conclusion 
is quite dependent on the assumption made for the maximum 
quantity of helium in the ejecta, that is still compatible with a 



type Ic event (above numbers are obtained for a maximum value 
ofO.5M ). 

Can we conclude that single-star models would produce 
most of the observed WR stars, while close binaries would pro- 
duce the majority of type Ibc? We think that this is likely too 
schematic a representation of reality. First, not all single-star 
models will produce a BH (even with the weak mass-loss rates 
used in the present work) and therefore some single stars will 
produce a visible type Ibc event. Second, in the absence of any 
reliable understanding of how massive stars explode (particu- 
larly if rotation is involved), it may be premature to draw firm 
conclusions regarding the nature of the visible counterpart of the 
final core-collapse. Third, the uncertainties concerning the mass- 
loss rates have a strong impact on both single and binary scenar- 
ios and still prevent definitive conclusions from both channels. 
Fourth, it would be striking that the predicted ratios for type Ibc 
SNe from single-star models would fall so well in the range of 
the observed values just by chance. Finally, studying the distri- 
bution of WR stars with respect to their host light distribution, 
iLeloudas et alj d2010h found similar trends between the distribu- 
tions of WN and WC stars compared to those of SNe lb and Ic, 
supporting the idea that WR stars are the progenitors of at least 
part of the type Ibc SNe. For all these reasons, we think that 
any strong statement attributing all Ibc events to close-binary 
evolution or to single stars is quite premature at present. Nature 
probably includes both channels, but with a frequency which still 
remains difficult to assess quantitatively. 

In this context, it is very interesting to mention two re- 
ce nt obser vations regarding progenitors of core collapse SNe: 
P lSmarttetail (12009) found no RSG type II-P progenitors in 
the mass range between about 18 and 25 Mq. In our opinion, 
this may well be explained by a strong mass loss during the 
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Fig. 10. Track of the rotating 20 Mq model in the B - V vs M v 
plane. The end point of the track is the blue point. The red square 
correspon ds to the estimated p osition of the progenitor star of SN 
2008ax (ICrockett et al.ll2008l) . 



RSG stage, which would lead stars in this mass range to ex- 
plode as a type II-L, type ITb, type lb, or type Ic SN (see 
Table 0}. Of course, the strong mass loss could be stimulated 
by the presence of a companion, but in that case one should pro- 
vide an explanation for why it only occurs above a certain mass 
limit. At the moment, we tend to favour an explanation based 
on some physical process becoming active o nly above a mass- 
lumin osity threshold, and related to pulsation ( Yo on & Cantielld 
2010) or supra-Eddington luminosity (Paper I). 2) Observations 
of yellow supergiant progenitors for lib SNe have been made 
recently (ICrockett et all 120081 see lMaundetaTll201ll and ref- 
erences therein). These progenitors likely evolved back from a 
RSG stage, and might be explained by some enhanced mass loss 
during the RSG stage. This supports the sta tement above and is 
well in line with the maser observations by Dav ies et al.l {2008) 
di scussed in Se ct. [4] This possibility has recently been discussed 
bv lGeorgvl(l2012l) . 

Let us mention here the particularly interesting case of SN 
2008ax discussed in Crockett et al. (2008). The authors mention 
the possibility that the progenitor of this SN could be a WNL 
star, with a very low hydrogen surface abundance. However, 
their non-rotating models do not reproduce the position of the 
observed progenitor star in the HRD well. Their best estimate is 
a ~ 27 - 28 Mq model. The difficulty of this model is the incon- 
gruously high mass of their CO core (around 9 Mq, while the 
light curve and ejecta velocities of the SN ar e compatible with a 
core mass of about 4-5 Mq; see lCrockett et a l. 2008). As shown 
in Fig. [TUl the track of our rotating 20 Mq model in the B - V vs 
My planqj ends noticeably close to the SN progenitor position. 
Moreover, it ends as a WNL star (see Table [T), with a H-content 
in the ejecta of 0.02 Mq, and a CO-core mass of 4.73 Mq (see 



5 The conversion of the quantitie s log(L) and lo g (r e ff/K ) to M v 
and B - V is done accordin g to [B oehm-Vitense ( 198 l|). | Flowed 
dl977h. ISek iguchi & Fukugltal d2000h . iMalagnini et all i 19861) and 
ISchmidt-KaleiU l982'). 



Table©. It is thus an excellent candidate for explaining the ob- 
served properties of the progenitor of SN 2008ax! 

6. Rotation rate of young pulsars and progenitors 
of LGRB 

Figure QT| shows how the specific angular momentum varies 
inside the rotating models at different evolutionary stages. It 
should be kept in mind that in the absence of any transport of an- 
gular momentum, the specific angular momentum would remain 
constant with time. In the current models, of course, transport 
mechanisms triggered by convection and rotation are present 
and therefore change the specific angular momentum. As a nu- 
merical example, we consider our 60 Mq model (right panel). 
At the beginning of the core H-burning phase, the specific an- 
gular momentum at the mass coordinate corresponding to the 
remnant mass is 3.2 • 10 17 cm 2 s . At the end of the core H- 
burning phase, the specific angular momentum at the same point 
is 3.6 ■ 10 16 cm 2 s , i.e. decreased by nearly a factor of 9. At 
the end of our computation, just after the core C-burning stage 
in this case, the specific angular momentum is 2.5 • 10 16 cm 2 s , 
i.e. decreased by nearly a factor of 13 with respect to the initial 
value. Since the evolution is very rapid up until the SN explo- 
sion, this number is likely to be a very good approximation of 
the specific angular momentum at the time of the pre-SN stage 
and at the border of the region that will be locked into the stel- 
lar remnant. Interestingly, it appears that more than 96% of the 
reduction of the specific angular momentum at that mass coor- 
dinate occurs during the MS phase. T his underlines, as already 
emphasised in previous works (see e.g. Hirschi et al. 2004, 2005; 
Woosley & Heger 2006), how important it is to properly account 
for angular momentum transfer during the MS phase. 

Another point worth noting is that a low surface velocity 
does not necessarily imply a low angular momentum for the cen- 
tral region of the star. Typically, at the end of the evolution of 
our 60 Mq model (see Fig. fTTT i. the surface equatorial velocity 
is only 22kms~', while the specific angular momentum at the 
edge of what will become the stellar remnant is well above the 
minimum required for a stable orbit around a Kerr BH, i.e. suffi- 
cient for allowing the accreting matter to form an accretion disk 
around the BH that will likely form at the end of the stellar life- 
time. 

Assuming that the angular momentum content in the mass 
coordinate that will become the remnant remains in the NS, we 
then obtain the pulsar period at birth indicated in Table 00. The 
periods for our models range between 0.08 and 0.13 ms. These 
are extremely short periods, not only much shorter than the ob- 
serv ed periods of young pulsars, which are between 20 and 100 
ms (iMuslimov & PageNl9"96l iMarshall et all 1 1998b iKaspi et alJ 
1994), but even shorter than the critical velocity of NS, which is 
between 0.44 and 0.65 m£] Some additional braking mechanism 
is therefore needed to reconcile our predictions with observa- 
tions (see below). 

A similar conclusion can be obtained by considering the 
frequency of LGRBs. Assuming that a LGRB is produced ev- 



6 To compute the periods, we assumed that all the angular momentum 
contained in the Lagrangian mass M rem , bar remains in the NS. We then 
used the relation I = 0.35 M rem . grav R 2 whe re I is the momentum of 
inertia of the NS JLattimer & Prakash 2001), and where a radius R = 
12 km is assumed for the NS. 

7 We estimated this critical period on the basis of the gravitational 
mass of our remnant, a radius of 12 km for the NS, and a critical 

GM 
R3 



(Keplerian) angular velocity given (in a classical way) by fl crit = _/ ^M- 
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Fig. 11. Specific angular momentum for the rotating 85 Mq (Left) and 60 Mq (Right) models as a function of the Lagrangian mass 
coordinate normalised to the actual total mass of the star. The continuous (black) curve corresponds to the ZAMS, the short-dashed 
(green) curve to the beginning of the core He-burning phase, the long-dashed (red) curve to the beginning of the C-burning phase, 
and the medium-dashed (purple) curve to the end of the C-burning phase. The dot-short-dashed (blue) curve indicates the minimum 
specific angular momentum for a stable circular orbit around a Schwarzschild black-hole and the dot-long-dashed (blue) indicates 
the minimum specific angular momentum for a stable circular orbit around a maximally rotating Kerr black-hole. The grey area 
cuts the lines at the point corresponding to the mass coordinate equal to the mass of the remnant indicated in Table [4] The surface 
equatorial velocities are indicated for each stage. 



ery time the three following conditions are fulfilled: 1) a BH is 
formed, 2) the angular momentum of the newborn BH is high 
enough to form an accretion disk, 3) a type Ic SN is associ- 
ated with th e event. These c onditions are those needed to obtain 
a collapsar dWooslevlll993l) . and to the observational associa- 
tion of broad-l ined Type Ic SNe with several nearby (z < 0.5) 
LGRB events dWooslev & Bloomll2006l: iBerger et al.11201 ll and 
references therein). With this hypothesis, we obtain favourable 
conditions for a collapsar in the mass range between 40 and 
60 Mq. In Table we show the expected ratio of collapsar to 
core-collapse SN (CCSN) obtained from present rotating mod- 
els, computed using Eq. dTJ. They are between 5 and 9% of 
the rate of core-collapse SNe. The fraction of LGRB to SN 
Ibc is much higher tha n the recent estimate of this fraction by 
Soderber g et al.l(l2010l) . which is about 1%, and is obtained from 
a radio survey of type Ibc SNqfJ 

From the above discussion, it appears that some additional 
braking mechanism is at work. What could this additional brak- 
ing mechanism be, and when does it intervene during the evo- 
lution of the star? Such a braking mechanism could intervene 
during the nuclear lifetime of the star, at the time of the core- 
collapse event and/or during the very first phases after the birth 
of the neutron star. 

At the moment, the solution that has been extensively ex- 
plored in the literature invokes a strong coupling mechanism 
active inside the star during its nuclear lifetime. This coupling 
could be caused by a strong interior magnetic field imposing a 
solid body rotation law during the core H-burning phase. The 
interior magnetic field could be fossil or generated by a shear- 



dynamo in the radiative layers. Coupled with the mass loss 
prescriptions used in stellar evolution, these models can repro- 
duce the shortest observed periods of young pulsars (Heger et al. 
2005) and restrict the progenitors of LGRB to those stars be- 
ginning their evolution with a very high initial velocity, so high 
that t hey follow a nearly chemical homogeneo us evolution (see 
e.g. lWooslev & Hegeril2006l: I Yoon et al.ll2006l) . Moreover, these 
models obtain favourable conditions for LGRB only at low 
metallicities. 

This solution is quite elegant in the sense that the physics of 
stars becoming pulsars (or NS) and stars becoming LGRB is the 
same (all have nearly solid body rotation during the MS phase). 
Only the initial velocities are different. However, there are also 
problems with this solution. First, there are still difficulties in ac- 
counting for young pulsars with the longest periods, indicating 
that even in that case an additional braking mechanism may be 
missing. Second, observations from asteroseismology indicate 
that a strong gradient of the angular velocity seems to be present 
in some massive stars (even stronger than in the present models), 
showing that at least some massive MS stars avoid solid body ro- 
tation ( Aerts 2008). Third, some authors (Zahn et al. 2007) have 
challenged the shear dynamo theory presently adopted in the 
models responsible for the strong coupling. Fourth, if the field is 
fossil, one can wonder why it has no surface components. A sur- 
face magnetic field may in the latter case produce a very efficient 
braking ( Mevnet et al.ll201 H) . which would make the formation 
of a collapsar nearly impossiblqj. 



8 This technique allows the detection of undetected LGRBs, since 
any relativistic outflow will produce a long-lived non-collimated after- 
glow. 



9 It would be extremely interesting in that respect to conduct spec- 
tropolarimetric observations to measure the geometry and the strength 
of the surface magnetic field, coupled with asteroseismological obser- 
vations giving information on the internal variation of the angular ve- 
locity. This would, for instance, answer questions such as whether the 
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Fig. 12. Evolutionary tracks in the HRD. The colours along the tracks indicate the minimum magnetic field (assumed to be aligned 
with the rotation axis and dipolar) able to couple the stellar wind with the stellar surface and thus to exert a torque. Higher magnetic 
fields than the minimum value indicated here woul d reduce the angular m omentum content of the star. The minimum value is 
obtained by imposing that the parameter 77 defined in Ud-Dou la et al.l (|2009) is equal to one. Left panel: Main-Sequence and RSG 
phases. Right panel: WR phase. 

Table 5. Mass range and ratio of LGRB compared to various SN types. 



mass range 



BH -> bright SN 
LGRB / CCSN LGRB / SN Ibc 



BH -> no bright SN 
LGRB / CCSN LGRB / SN Ibc 



40.0 - 60.0 M e 
36.0- 67.0 M e 



0.049 
0.076 



0.254 
0.394 



0.056 
0.087 



0.663 
1.029 



From the points listed above, we are inclined to think that 
most of the massive stars do not rotate as solid bodies during 
their MS phase but possess a differentially rotating interior. If 
true, this implies that the slowing-down of the core needed to 
reproduce the observed periods of young pulsars should result 
from another process than the strong coupling caused by an 
internal magnetic field. Could it be caused by a surface mag- 
netic braking? A sufficiently strong surface magnetic field may 
indeed force the wind material to follow the magnetic lines. 
Conservation of an gular momentum will exert a torque and thus 
slow down the star (Ud-Doula et al . 2009O It does appear now 
that at least some massive stars show a surface magnetic field 
(see the review by Petit 2011). We show in Fig. [12] the mini- 
mum value of a dipole-aligned surface magnetic field needed to 
interact with the stellar winds. We see that on the MS, for stars 
with masses between 20 and 40 Mq, even a < 200 G magnetic 
field would be sufficient to impose some coupling between the 
wind and the surface and would therefore slow the star down. 
The minimum magnetic field allowing a coupling with the wind 
is even lower during the RSG phase, with a couple of 10 G. In 
stars more massive than about 50 Mq, magnetic fields between 



stars presenting a strong surface magnetic field rotate as solid body in 
the interior. 

10 The magnetic field could also drive the matter back to the surface 
of the star and thus prevent any loss of angular momentum. In that case, 
the magnetic field would not help as a braking mechanism. 



300 and 600 G are needed during the MS phase, and are of the 
order of a few 10000 G during the WR phase. The consequences 
of such a surface magnetic field braking mechan ism remain to 
be studied (see first results in Meynet e t al.ll201 lb . 



We finally note that a large difference between models with 
solid-body rotation (and thus "strong" internal coupling) and 
those with internal differential rotation ("moderate" internal cou- 
pling), such as those presented here, is the manner in which 
LGRB occurrence may depend on metallicity. For strong cou- 
pling, the occurrence of LGRB is mostly restricted to low metal- 
licity where the mass loss by stellar winds is weak. For moder- 
ate coupling a metallicity dependence is also present, but LGRB 
could occur at solar metallicity, as illustrated by t he present mod- 
els, and probably at higher than solar metallicity (Lev esque et al.l 
2010a b). Very interestingly, the lack of correlation between the 
energy of the explosion o f LGRB with metallicity as obtained 
by iLevesque et al.l (l2010ch may reflect the fact that the internal 
conditions, which govern the energy, are not strongly impacted 
by metallicity. This somehow supports the "moderate" coupling 
scenario, in which the properties of the core, in terms of angular 
momentum content, are much less sensitive to metallicity than 
in the "strong" coupling scenario. 
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7. Synthesis and conclusion 

We have presented predictions of single-star models and com- 
pared them to several important observational features: the pop- 
ulations of WR stars, of type Ibc SNe, and of collapsars at solar 
metallicity. 

We confirm with the new generation of rotating models 
and the results obtained from pr evious work (tMevnet & Maeder 
120031 120051; iGeorgv et al.ll2009l) that rotation favours the forma- 
tion of WR and type Ibc SNe. 

On the basis of these new models, we estimate that the 
single-star channel could account for about 60% of the observed 
WR stars at solar metallicity, with the remaining 40% likely orig- 
inating from close-binary systems. 

The fraction of type Ibc SNe coming from the single-star 
channel will remain difficult to obtain in a reliable way until SN 
models are able to reliably predict the properties of the explosion 
of core-collapse events leading to BH formation. At the moment, 
all we can say is that at least ~50% of type Ibc SNe could origi- 
nate from the single-star scenario. 

The present rotating models predict a lack of type II-P pro- 
genitors above ~ 17 Mq while the non-rotating ones predict 
such a lack only above 19 M<7). Both mod els therefore agree 
with the observations bv lSmartt et al.l (2009). who observed no 
type II-P SN progenitors with initial masses above about 18 Mq. 
Moreover, rotating models provide an excellent explanation for 
the characteristics of the progenitor of the SN 2008ax. 

We also discuss some possible shortcomings of the present 
models. Part of these shortcomings may also be due to uncer- 
tainties in the observations, or to a presently biased view linked 
to incomplete sampling. Based on the available observations, it 
does appear that the current rotating models may overestimate 
the duration of the WNL phase and underestimate those of the 
WNE and WC phase. A possible reason for this is the fact that 
the pre sent recipe for the WNL ma ss-loss rate predicts very low 
values dGrafener & Hamannll2008l) . 

Concerning the origin of the low-luminosity WC stars, the 
results may be in part improved by an enhancement of the mass- 
loss rates during the WNL phase. An enhancement of the mass- 
loss rates during the RSG phase and/or mass transfer in close- 
binary systems could also help resolve this problem. 

The question of the rotation rate of pulsars indicates some 
additional angular-momentum loss from the central region with 
respect to the scenario obtained in our models. While a strong 
internal coupling (caused, for instance, by an internal magnetic 
field) might resolve this problem, we mentioned alternative so- 
lutions. In particular, the role played by magnetic braking at the 
surface needs to be explored in more detail in this context. 

From a theoretical point of view, the conditions for hav- 
ing collapsars are not necessarily restricted to low metallicity. 
The effect of metallicity on the evolution leading to a LGRB is 
strongly correlated to the degree of coupling between the core 
and the envelope in stellar models. A strong coupling is more 
prone to confine GRB occurrence to low metallicity than a mod- 
erate or a weak coupling. 
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